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AN AUTOMATIC MACHINERY FAULT DIAGNOSTIC 
METHOD AND APPARATUS 

FIELD OF THE INVENTION 

The invention relates to an automatic fault diagnostic procedure for machines, 
such as steam, gas and hydro-turbines, generators, pumps, electrical motors, etc. 

BACKGROUND OF THE INVENTION 

Various attempts to monitor the health and make fault diagnoses of machines 
have been attempted. Conventionally, this is done by human experts after abuoimal 
behaviour, such as increased vibration level, has been observed- The experts may in 
cases be supported by expert systems. 

The fault diagnosis may include analysing vibration and process data jfrom 
critical machines to determine whether the specific machine components show a 
tendency of unbalance, misalignment, cracks^ wear or other faults. 

As regards the use of expert systems, the results &om such systems may 
merely indicate the association of a symptom with the likelihood of a potential fault, 
but does not predict the future development of this fault. Fiutheraiore, expert systems 
are sensitive to boundary conditions of the collected data and knowledge, and may 
even indicate faxilts, which would possibly not occur in practice. In addition, such a 
method, when monitoring the many machine components, soon arrives at a point 
where the amount of data to be analysed reach the limit of the capacity of the 
computer executing the algorithms of the expert system. 

Expert systems and other types of automatic fault detection systems also 
require training, and the difficulty of training an expert system for the many different 
types of fault and indicators thereof may be prohibitive. 

There accordingly remains a need for an improved automatic machinery fault 
diagnostic metliod and apparatus. 

SUMMARY OF INVENTION 

According to a first aspect of the invention there is provided: a method of 
automatic fault diagnosis for machinery having a plurality of componi^ts, based on a 




comprismg: for at least one machine component, and for at least one fault that may 
occur in that component, cairying out the steps of: calculating for each of a plurality 
of symptoms for indicating that fault, a symptom value as a function of the reduced 
dataset calculated from measured data; combining the symptom values to give a total 
fault symptom strength value indicating the strength of the fault; and carrying out 
faiilt processing if the total fault symptom strength lies above a predetermined value. 

By calculating parameters representing the strength of a number of physical 
symptoms, and then combining these parameters to give a value of total fault 
symptom strength^ it is possible to automatically evaluate the state of a component 
with regard to each of a nmnber of possible faults. This approach works much better 
than expert systems or approaches using automatic training of a system using 
traiiaing data because the method according to the invention is able to focus on the 
specific featm-es of the dataset indicative of each fault. 

Furthermore, by talcing into accoimt a classification of the machine state the 
method may only use data from relevant machine states. For example, data taken 
when a machine is idling may be irrelevant to the case when the machine is operating 
at full power. 

The stqp of combining the symptom values may include calculating the sum 
of the symptom values raised to the power of a predetermined constant and taking 
the predetermined constant root of the sum. In this way, it is avoided that a single 
outlier value has too great an effect on the final calculation, while still combining all 
the symptom values. 

Tht step of combining the symptom values may include multiplying the 
combined symptom values wifli a value representmg the strength of tiie symptoms. 

The signature dataset may be a database. The method may include 
automatically taking measurements of the component when predetermined 
conditions occur, such as at regular intervals, and testing whether the new 
measurements r^resent a change on the reduced data stored in the signature 
database. If there is no change;, preferably simply this fact is recorded, alternatively 
the data may simply be abandoned* This approach avoids adding uimecessary data to 
the signature database tliereby reducing the amount of data required to be stored. 

The fault processing may include executing different program code 
depending on the calculated total fault symptom strengtli. 
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For example, if the total fault symptom strength lies below a first 
predetermined value, no further processing may take place for that fault. This allows 
fault symptom values which nxight have arisen by chance to be ignored. 

If the total fault symptom strength lies above the first predetermined value the 
5 method may, with respect to that fault, switch from a first mode in which the value of 
the total fault symptom strength is not recorded on an ongoing basis to a second 
mode in which the total fault symptom strength is recorded on an ongoing time series 
basis. The recordal of the value of the total fault symptom strength each time tlie 
fault system strength is calculated allows the time evolution of the fault to be 
10 recorded. 

The fault processing may further include outputting suitable warning 
messages, for example on the first occasion that a fault has been detected, and/or 
when the total fault symptom strength exceeds a second predetermined value 
indicating that the fault is seiious or critical. 

1 5 The method may include predicting the lilcely value of tlie total fault 

symptom strength at vaiious fixture times, based on the time series. Warning 
messages may be output when the time series indicates that the total fault symptom 
strength will exceed a predeteimined value after a predetermined time. 

The method may include predicting the fioture values of Ihe total fault 

20 symptom strength by automatically fitting the time series data to a plurality of 
different models and using the best fit. For example, the time series may be 
modelled as a straight line or as a polynomial of given order and the best fit taken. 

The method may include automatically calculating the faults for each 
component when predeteraiined conditions occur. Conveniently, the predetermined 

25 condition may include the elapse of a predetennined time after the calculation was 
last performed whereby the system monitors the health of the components on a 
regular basis. 

Preferably, the method includes measuring machine data and updating the 
reduced database on an ongoing basis. 
30 In another aspect, the invention relates to a method of automatic fisiult 

diagnosis including for each of a plwality of components of machines^ carrying out 
the steps of; for at least one possible fault in the component: calculating a total fault 
symptom strength value fi:om measured d&ta indicating the strengJii of the fault; and. 



if toe total fault symptom strength value lies above a predetermined value for that 
fault, recalling stored data relating to the total fault system strength value as a 
function of time, fitting the stored data to a trend line and predicting the time when 
the total fault symptom strength value will exceed a predetermined value, selecting a 
message based on the total fault symptom strength value and the predicted time» aud 
outputting the selected message. 

In this way the time to failure of the component may be estimated and useful 
messages provided to an operator to allow maintenance to be planned. 

The invention also relates to a computer program system for carrying out the 
method. The computer program may taclude a number of program blocks recorded 
on a data carrier, the program blocks being arranged to run on various computers of a 
computer system including a number of computers, or on a single computer. 

The computer program system is preferably implemented in an object 
oriented language. A machine component object is preferably provided for each of 
tlae machine components for which autodiagnosis is performed. Preferably, the 
machine component objects are specific instances of a general object class for the 
type of machine component. 

This approach brings real benefits in setting up the system for a particular 
plant. A computer system may contain general classes of object, such as a shaft, a 
bearing and a foundation. Each of the compon^ts of a piece of machinery can be a 
child of the relevant general class, inheriting the general properties of the class but 
furflier programmed with the specific features of that partioular component. 

In particular, the general class of object may include members for each of tlie 
different possible faults for objects of the class. 

The invention also relates to a system for automatically outputting a fault 
diagnosis for at least one possible fault, including code for carrying out the method 
discussed above. 

In machines used in power plants, there are often several hundred critical 
components, and since several machines are in use simultaneously, and these 
machines further have to be keep imder surveillaace snnultaneouslyj the amount of 
data from the involved measuring transducers increases very quickly. Embodiments 
of the invention achieves a reduction of the observed data. 

Furfhennore, ttie invention may take into account particular behaviour of the 



individual critical machine components during different operating conditions prior to 
analysing acquired data therefrom. In this way, the measured data may be 
interpreted more completely. Since different specific features are required from the 
acquired data, the "raw" data may be repeatedly subjected to altemative signal 
processing procedures to enhance particular fault features before the results are 
interpreted. The measured data may, for instance, be subjected to windowing, Fourier 
tt'ansformation, statistical estimation, filtering, integration, differentiation, wavelet- 
transform or alternative signal processing methods, resulting in a fault specific 
reduced dataset or "signature". The reduced dataset may then be interpreted, using a 
particular Method. The signature is the result of a sequence of signal processing 
functions that transforms the observed data to a domain that enhances fault specific 
features. Specific obsei-ved parameter values, such as temperature or vibx^ation level 
may be insignificant in one operating state and, in another operating state, may 
indicate a critical condition for the machine component . The actual operating state 
of the machine should accordingly be considered, before the measured data are 
interpreted. 

The procedure, according to the invention^ may use a machinery fault class 
library including references to specific signatures, i.e. specific pieces of data in the 
dataset, calculated from signals acquired from sensors placed at specific locations on 
said machine. 

The analysis may be performed on the heavily reduced data sets being stored 
in the server database. As a result, the time it takes the auto-diagnostic software to 
process the stored data is reduced considerably. Furthermore, since the auto- 
diagnostic software in addition works with signatures specific to particular kinds of 
faults in particular components, lixe risk of misinterpretations of the results and, 
hence, the untimely issue of false alarms are also reduced considerably. 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understandmg of the invention, embodiments of the invention 
will now be described, purely by way of example, with referesnce to the 
accompanying drawings in which; 

Figure 1 illustrates a system according to the invention; 

Figure 2 is a flow diagram illustrating processing for a single fault; 



Figure 3 is a schematic diagram of fault processiag overall; 
Figure 4 is an inlieritance diagram for component classes; 
Figure 5 is a schematic diagram of the configuration Wizard; 
Figure 6 is a schematic flow diagram relating to updating data- 
Figures 7 to 1 1 illustrate a first symptom of rotor unbalance; 
Figure 12 illustrates output of the configuration Wizard; and 
Figure 13 illustrates the structure of the Faxilt Library. 
Like or correspondiug components are given the same reference numerals in 
the different figures* 

DETAILED DESCRIPTION 

Referring to Figure 1, a system according to an example of the invention 
includes a computer system 2 having a network 4 connecting a plurality of 
workstations 6, a server 8 and a signal processing unit 10, 

The signal processing unit 10 is connected to a number of sensors 12 on a 
machine component 14, here a rotor. In practice^ sensors 12 mil be provided on a 
number of components 14, optionally on a number of machines 16 which may have a 
plurality of components 14. A number of sign^ processing units 10 may be provided 
if convenient. 

The server 8 includes at least one processing imit 18 and a memory 20, the 
latter storing program code 22 for automatically diagnosing system faults. . The 
sldlled person will appreciate that the memory may include RAM, EEPROM, hard 
and floppy drives, DVD and olher data storage means. 

Program code may also be provided on workstations, signal processing unit 
10 or indeed on any other computer connected to the network. 

A Fault Diagnostic System of program code 22 runs on the central server 8 
and optionally also any mmiber of clients 6 (operator/user temiinals), which are 
connected through local or global networlcs 4. The central server 8 uses a general 
database management system for sorting, classifying and processing information or 
for activating an alarm when a fault occurs* The code 22 uses a particularly 
developed scientijSc diagnosis language syntax (see appendix A) and the method 
fomialises the fault features as seen in the fault selective signatures. The semantics is 
founded on physical mathematical models and empirical knowledge for the particular 



faults and components. 

During operation, the central database server 8 executes simultaneously 
several processes such as a data server, a user interface, an alarm server, graphical 
data presentation, etc. Most important for tiie present patent application is an 
5 " AutoDiagnostics Inference Engine'* (auto diagnosis server), a process canning out 
the diagnostic method for all components of all machines connected to the Fault 
Diagnostic System.. 

The or each signal processing unit 10 measures data from the sensors and 
updates data stored in a measurement database 44 on the central server on a regular 
10 basis by monitoring inputs on the sensors to obtain measured sensor values, as 

illustrated in Figure 6. Measurements are taken periodically (step 60). Fault specific 
reduced data is calculated from the measurements (step 62) for each of the possible 
faults , The reduced data is then compared (step 64) with data stored in 
measurements database 44, and if there is a significant change the data in database 44 
15 is updated (step 66), Otherwise, the new reduced data is discarded, hi this way, vast 
amounts of identical data are not stored in the measurements database 44, thereby 
keeping measurements database 44 manageable. 

The measurements are used to determine the classification of the machine 
state and the reduced data is compared with data stored for the same machine state. 
20 The machine states may include idlyingj run-up, full power^ or other states as 
appropriate. The reduced data is stored in the database 44 together with the 
classification of the machine state. 

The steps carried out by the autodiagnosis code 22 in the central server 8 will 
now be described with reference to Figures 2 to 5, 
25 An "AutoDiagnosis Supervisor" process runs transparently in the background 

of the system and causes fault calculation to occur for each component and each fault 
and issue messages to the user, as illustrated in Figure 3. As shown in the pseudo 
codcj the AutoDiagnosis Supervisor process runs througli all machines, all 
components and all faults to get individual methods firom the fault library, 
30 component parameters and signatures firom the configuration database, and executes 
the diagnosis in accordance with the method. The AutoDiagnosis Supervisor process 
may be initiated either on a fixed schedule or when triggered by an event, for 
example an alarm. In the case that faults are detected messages and 
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recotmnendations are gathered from the databases and issued to all clients (users) on 
the system. 

The processing for each fault may be run sequentially for each machine, each 
component, and each fault, as indicated schematically in Figure 3, Alternatively, it 
5 may be preferable to commence a separate process for each component, or even each 
feult, when predetermined conditions occur, such as at regular intervals. In this way, 
some components can be checked more frequently than other components- 
Figure 3 also illustrates the storage used at various pomts of the metiiod. 
Reduced dataaets obtained from measurements are stored in database 44. A further 
10 database 46 stores flie configuration of the equipment, and a yet further data store 48 
stores a fault library, including the code for each possible fault. All tiiese may be 
stored in memory 20 of the server 8 or elsewhere as convenient, 

Figure 2 illustrates the calculations carried out by code 22 for a specific fault 
on a specific component. The specific fault may include unbalance of a rotor, 
15 misalignment, instability, and various other possible faults that will depend on the 
specific component. 

The code for the specific fault is started in step 24, Then, in step 26 symptom 
values for each N possible physical symptoms of the fault is calculated, where N is a 
positive integer. 

20 An example of how this works is presented below for the case of one possible 

fault, unbalance of the rotor. In this example, the symptom values in this example 
include symptom strength values v,- which indicate the relative strength on a scale of 
0 to 1 of the deviation firom a baseline value measured in the component when 
operating the same class of operation and symptom severity values Wi which are 

25 absolute values. 

Depending on the fault and the symptom, data is selected from the 
appropriate classification or classification of machine states, Thus^ for some fault 
symptoms, data from one machine state will be relevant, whilst other fault symptoms 
may require data from another machine state. The symptoms take data from the 

30 signatures (the reduced data stored in database 44) from one or more classifications 
of machine state data as appropriate. Each symptom cycles in turn between data 
fi-om different machine states if required. 

Next, in step 28, the symptom values vi. are combined to give a total fault 
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symptom strength value T. Iq the specific example, T may be given by 
T - Fuzzymin(wi ... w„) . Fuzzymax (vi ...vO 

5 The flanction Fuzzymin gives a fuzzy ijiinimw of its arguments and 

Fuzzymax gives a fuzzy maximum of its arguments. 

In step 30, it is determined wheflier to take any action. If the total fault 
symptom strength value is below a predetermined value, for example in the range 0.3 
to 0.4, no fiirtlier action is taken with regard to that possible fault, and the processing 
10 for the fault ends (step 32). This allows small, probably spurious fault values to be 
ignored. 

The total fault symptom strength value history is tlien recalled (step 34), and 
a prediction made based on the history as to when the fauh becomes serious (step 
36). Depending on the total fault symptom strength value and the prediction, one of 
fl 15 a possible number of messages is selected (step 3 8) and sent (step 40). The new total 

fIJ feult symptom strength value is added to the symptom history (step 42). 

A message may be sent, for example, the first time the fault is detected, and 
then again when sufficient history data has been collected to malce a prediction with 
confidence limits less than a predetennined value. An alarm message may be sent 
O 20 when the total feult symptom strength value exceeds a predetermmed value. 

The method provides a total fault symptom value, a scalar value in the range 
O-l^ and is classified in three degrees of severity; Fault Detected, Fault Prediction 
and Critical Fault. Take as an example, the fault of "wear and tear", becoming 
detectable, Le. exceeding the sigmficance level of detection. At an early stage, the 
25 fault development may be far j&om critical and the fault is simply detected, without 
tite system having enough symptom data to provide an unambiguous prediction of 
the fault development trend. After some time, the fault may tend to develop in a 
significant yet predictable manner, as function of time. The system is then capable 
of predicting how flie fault develops in the future, and, furthermore, estimating the 
30 date at which the fault will be critical. As die time for predicted and critical faults 
become available, the operator has the opportunity to performing scheduled plant 
maintenance/corrective activities and enjoy the benefits of preventing nnscheduled 
plant outages. A critical fault is a fault that has developed to a level considered 
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dangerous for continued operation. The occurrence of such faults nearly always 
implies an immediate closedown of the machine in question for repair or exchange of 
defective components. Any of the three fault levels cause alarm messages when first 
detected with an indication of the machine^ the component, the type of the action to 
be takexa together with the prediction of future fault symptom developments, if 
predictable. 

An important feature of the system iis that it provides output messages clearly 
identifying the faxilt, taken fiom the message database in the data store 48. In this 
way, the operator can clearly identify the fault that the message relates to. 

A plant may have several machines, each machine m^iy have several 
components^ each component may have several potential fiaults, each fault may have 
featm-es to be seen at many signatures, but each fault/signature component lias only 
one method to determine the fault symptom. The AutoDiagnosis supervisor 
(inference engine) gathers for each faulty measured signatures fi-om the measurement 
database^ method firom the fault library and executes the method to deteimiue partial 
symptom strength. 

The fault library 48 is organised in a hierarchy that reflects the plant 
configuration as illustrated in Figure 13. It will be appreciated that in the object 
oriented approach used in the examples the objects stored in the fault library may 
include the code needed for calculating and detecting the corresponding fault. 
Code for causing processing of each of the varioiis feults are stored in the fault 
library 48, as objects relating to the different components, the faxilts for each 
component, the signatures for each fault and the methods for each signature. 

The fault detection code 22 stored in the feult library 48 is conveniently 
programmed in an object oriented programming language with inlieritance, A 
separate computer program object may be provided for each machine component. 

Figure 4 illustrates the inheritance of these separate computer program 
objects 52. A plurality of general fault class objects 54 are provided, each 
corresponding to a type of component such as a rotor, a bearing, or any other 
component. The possible faults for that type of object are encoded in the general 
fault class objects 54. Each specific object 52 inherits from one of these general fault 
class objects, but further includes component specific information such as size^ 
intended running speed and tbe like, 
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Conveniently, liie separate computer program objects 52 further inlierit from 
a class for each machine 56 and each plant 58, 

By programmmg in this way it becomes very easy to update the system to 
detect faults in new objects, Ahnost all of the code is provided in the general class, 
and this code needs little alteration for the new component. 

Indeed, the new code may be generated by a wizard as illustrated in Figure 5. 
The AutoDiagnosis Fault Library Wiz;ard provides tools for building and editing the 
underlying Fault Class Methods. Hie information is accommodated in data store 48 
containing all information regarding the machine con^onent, component specific 
parameters, fault description and alarm messages, descriptive information on how the 
signatures shall be measured and the diagnostic class method. As the pseudo code in 
the figure suggests, the wizard requires user input ("Answers'*) relating to features of 
the component ("Questions") and the user input is used to update the configuration 
database 46 for every machine and component in the plant, An example user input 
screen is shown in Figure 12. 

For every potential component fault, specific measurements, as already 
configured for the machines, are associated with this fault and, as impUed, with a 
fault class method in tlie fault library. 

For a better understanding, a glossary of terms used herein follows. 

Fault: A fault is a mechanical deviation of the component, that endanger the 
intended functionality, and which has, may have or will have an impact on the 
reliability or efficiency of the component, if it develops further- 
Signature: A signature^ is the result of a sequence of signal processing 
operations that transforms the observed data into a domain that enhances fault 
specific features, such as autospectrum and the Selective Envelope Spectrum^ which 
reveals changes in ftequency components, and detects amplitude modulations and * 
impulses in a particular frequency range respectively. 

Fault Class Method: The Fault Class Method is a formalised diagnostic 
method to detemiine the Symptom Strength for a particular fault on basis on one or 
several signatures acquired from a particular machine and machine component. The 
Method itself is a procedural diagnostic syntax of a scientific diagnostic language, 
developed specifically to provide a rational and focused semantics for determining 
faults from acquired signatures. Syntactic elements of an example Fault Class 
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Method ai'e listed in Appendix A. 

Symptom Strength; reflects the degree of development state of the individual 
faults as seen from the associated signatures* 

System Elements: The system has a set of application elements that provides 
5 functionality for establislmient and maintenance of the library of Fault Class 

Methods^ to configure Fanlt Class objects, viz. a representation of a physical plant of 
machines/components/faults and associated measurements and a process that 
executes the Methods on the acquired signatures: 

The system may typically be installed and operated on a plant in according to 
10 the following procedurej Initially, the critical machines 16 and their constituent 
components 14 to be kept under surveillance are selected. Subsequently, a plurality 
of transducers or sensors 12 are mounted on the components 14, for instance 
accelerometers and proximity probes (for absolute and relative vibration 
81 measurements respectively), tacho probes (for measuring rotational speed), 

1 5 thermocouples (for measuring temperatures) flow sensors (for measuring liquid 
f y flows), gap probes (for measuring relative distance), etc. These transducers 12 are 

connected to one or several signal processing units 10 (SPUs), capable of 
CI quantisation, processing and caching the information from the individual transducers 

east . 

[ and eventually provide data commimication to the central system server via 
Q 20 standard networks, A typical plant, such as an electric power plant or a cement 

manufacturing plant may have several himdred transducers 12. 

EXAMPLE 

The principles of conJiguratLon and execution of auto diagnosis will now be 
25 illustrated for a typical machine component; a rotor. In the following, one practical 
example of configuration and operation of the Fault Diagnostic System will be 
e>cplained- The example concerns a rotor supported in journal bearing. The system's 
Machine Fault Library includes three g^eral classes of rotors characterised by 
mechanical configuration and applicable instrumentation; one of these, the rotor with 
30 journal bearings, is selected 

Fig 1 shows a rotor 14 as part of a generic machine supported by two journal 
bearings 80. At each bearing 80, two displacement probes82 and one aocelerometer 
84 are mounted. The signals &om the displacement probes 82 and the accelerometer 
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84 are fed to a signal processing unit 10 for calculating the particular signatures, that 
subsequently ai'e send to the central server: 8 and saved in the database 44. 

The potential faults for the rotor are unbalanoej naisalignment, bowed shaft, 
cracked shaft, bearing instability, mechanical looseness and rotor rub. 

In the example, consider the case of one of these faults " unbalance". 4 
signatures per 6 sensors (totally 24 signatures) are used to determine five partial 
symptoms^ here distortion, local phase, roundness, phase trend and to^-law. These 
partial S3miptotQS are combined in a flizzy logic sense and scaled by the severity of 
the vibration level to provide tlie current Total Symptom Strength for ttxe fault; 
unbalance. 

The three transducers 82,84 on each of the two bearings S6 supporting the 
rotor 14 comprise, as previously mentioned, two radial mounted displacement probes 
82, measuring relative displacement between shaft and bearing house and one 
accelerometer 84 measuring the absolute acceleration of the bearing house, If the 
displacement probes are moimted at an angle which to one another, with respect to 
the axis of the rotor, which differs from 90*^, the system ensures that the signals from 
said probes are geometrically transformed into signals corresponding to an angle 
separation of 90°. Another displacement probe provides a rotor angular reference 
signal for Order Tracking Analysis and Shaft Speed detection. 

The sensor signals can thus be summarized as follows: 
Displacement probe, X direction^ Driving End, XDE. 
Displacement probe, Y direction, Driving End, YDE, 
Accelerometer on Driving End bearing, ACCDE 
Displacement probe, X direction, Non Driving End, XNDE, 
Displacement probe, Y direction, Non Driving End, YNDE. 
Accelerometer on Non Driving End bearing, ACCNDB 

MEASUREMENT STRATEGY 

All the symptoms for unbalance on rotors with journal bearings are calctilated from 
OTA, Order Track Analysis, The frequency range should be selected to allow OTA 
to produce measurement for the first four multiples of the rotating frequency* OTA 
should be recorded during nm-up, and OTA time history should be kept. 
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SIGNATURE 

All the symptoms for unbalance on rotors with journal bearings are calculated firom a 
signatm e known as OTA, Order Track Analysis. The amplitude$ of vibration 
measured by the sensors 82,84 for the first fonr multiples of the rotating frequency 
and the first order phase are used. Four kinds of information are used; current values, 
baseline, time history data and bode data ft'om last run-up. 

Two symptom values, which may be called a symptom severity and a symptom 
strength may be calculated for each symptom in the example. In other arrangements, 
only a single symptom value might be calculated. 

There are six symptoms on unbalance for rotors with journal bearings implemented: 
Symptom 1 . Dominating first order 

Symptom 2. Same vertical and horizontal phase shift over rotor 
Symptom 3. 90 degrees phase shift vertical/horizontal 
Symptom 4. Orbit roxmdness <5 
Symptom 5, First order phase steady versus time 
Symptom 6. First order proportional, to RPM^ during mn-up 

SYMPTOM STRENGTH 

The symptom strength for each of the symptoms of the fault unbalance is calculated 
by comparing the cuirent first order amphtude compared to baseline. As amphtude, 
the mis sum from all six setisors is used. The aocelerometer readings are integrated to 
displacement with use of the RPM before tlie rms summation- 
FaultStrength - loglO{(curr«it first order rms) / (baseline first order rms)} . 

The fault strengfli may also be considered to be a symptom strength for each of the 
symptoms of the feult, In embodiments, dififeent symptom strengths may be 
calculated for each symptom. 



SYMPTOM VALUES 

SYMPTOM 1. DOMINATING FIRST ORDER 

On the current OTA data, distortion is calculated. Distortion is defined as: 
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i^rx total rms of orders 2^3 

distortion = 100 — ■ — ; ; , ^ ^ ' 

total rms of orders 1, 2,3 

The rms is taken as a (squared) sum over all available seasors. The aocelerometer 

values at*e integrated into dtsplaoejnent with the use of the current RPM before the 

summation. A dominating first order will set the distortion to near zero, A smooth 

5 function fiizzygauss then transfers the distortion value into a symptom value between 

0 and 1. See figure I for the relation between symptom 1 and distortion. The 

symptom limit (symptom - 0.5) is initiaJly set to a distortion of 15%. 

SYMPTOM 2. SAME FIRST ORDER PHASE SHIFT OVER ROTOR FOR 
10 VERTICAL AND HORIZONTAL 
J DIRECTIONS, RESPECTIVELY 

From the current OTA data, the first order phase shift over the rotor is calculated in 
vertical and horizontal direction, respectively. A smooth fiizzygauss fimction is then 
applied on the difference between the phase shifts. The limit (symptom = 0.5) is 
^ 1 5 initially set to 0±30 degrees. See figure 2. 
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CI SYMPTOM 3, 90 DEGREES FIRST ORDER PHASE SHIFT BETWEEN 

VERTICAL AND HORiZONTAL 
DIRECTIONS 

20 From the cuirent OTA data, the first order phase shifts between vertical and 

horizontal directions are calculated for both sides of the rotor respectively. Then a 
vector addition, using die first order amplitude, is performed on the two sides to 
produce a ''total" phase shift. A smooth fixzzygauss function is then applied on the 
resulting phase. See figure 3. The limit (symptom = 0.5) is initially set to 90±20 
25 degrees. 

SYMPTOM 4, ORBIT ROUNDNESS <5 

From the current OTA data the first order orbit roundness is calculated. The major 
and minor axes of the orbit ellipse are calculated fi'om the first order magnitude and 
phase. Roundness is then defined as the ratio between the major and minor axes. Tliis 
30 means that a circle has roundness one, while a fl;attened ellipse has a greater 

roundness. (The term roundness is misleading, but is used in this feshion by many 
authors). The roundness is calculated for both bearings and then combined in a 
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smooth fashion by flie function teymax. Finally, a smooth function fuzzygauss is 
applied to get the symptom. See figure 4, Initially, the symptom limit (symptom = 
0.5) is set to a roundness of 5. 

5 SYMPTOM 5. FIRST ORDER PHASE TIME HISTORY STANDARD 
DEVIATION AND DRIFT 
< 10 DEGREES 

For all six sensors, the first order phase tune histories are loaded with GetHist. 
Initially, the histories ftom the last 20 days are loaded. The corresponding time 

1 0 vectors, which may be different for each sensor, are also loaded. The phases are 
unwrapped to avoid steps in the time histories, For a first partial symptom the 
standard deviations for all six phase time histories are calculated. The six standard 
deviations are combined with a smooth fu2Kymax function. The result is converted to 
the first partial symptom with a smooth fijzzygauss fimction; initial limit (symptom = 

15 0,5) is set to 10 degrees. 

Next, a linear regression is perfomied on the time vector and phase time 
history for each sensor. A fimction linreg calculates the linear regression coefficient, 
weighted with the square of the correlation coefficient. The coefficients for the six 
sensors are smoothly combined with a fiizzymax function. Then the total drift during 

20 the last 20 days is calculated and converted to the second partial symptom with a 
smooth fiizzygauss function; initial limit (symptom " 0,5) is set to 10 degrees. 
Finally, the two partial symptoms (phase standard deviation and phase drift) are 
combined in a smooth fuzzymin function. See figure 5. 

25 SYMPTOM 6, VIBRATION AMPLITUDE PROPORTIONAL TO RPM^2 FOR 
LOW RPM IN LAST RUN-UP 

For all six sensors, the first order amplitude versus RPM and the corresponding RPM 
parameters are loaded with the function GetBode, From the RPM parameters 
(RPMmiiL, RPMmax, number of RPMs), RPM vectors are calculated for each sensor. 
30 Initially, only RPM values up to 70% of the first critical RPM are considered* For 
those RPM values, a "Squareness Coefficient" is calculated with the function 
squareness. This coefScient tells how large part of the first order amplitude RPM 
history that is proportional to RPM^. The aocelerometer readings are integrated to 



17 

displacement before the Squareness calculation. The Squareness coefficients for all 
six sensoi-s are combined with a smooth Fu^zymin function. 

TOTAL FAULT SYMPTOM STRENGTH 

Finally, the Total Fault Symptom IStTength can be determined: 

TotalSymptomStrength = Fiizzymin({Symptom Value(l):(6)}) * 

Fnzsymax ({SymptomStrength)} 

Thus, if several of tlae symptoms are not present, the corresponding 
SymptomValues will be low leading to a low TotalSymptomStrength, Only if the 
symptoms are present will a high TotalSymptomStrength be obtained^ even in the 
presence of significant vibration. This makes it easier to for the method to 
distinguish the fault. 

This new parameter (TotalSymptomStrength) is saved as function of time and 
Machine State and is used in the Fault Strength Model for prediction of fault 
development rate and residual safe operating life for the component 

The skilled person will reaUse that the system is not limited to the example 
presented above. The system is capable of detecting, not just mbalance symptoms in 
a rotor, but also a variety of symptoms of a variety of types of fault in a variety of 
components in a variety of different types of plant 

In particular, the details of the program steps and flow charts may be 
modified whilst remaining withm the scope of the invention. 
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APPENDIX; ImplemOTted Method Functions 

To definine the Fault Class MeHiods, the fonnal sj^ntax of tlie Scientific Diagnostic 
Language provide tlie functions shown in the following table. 



binary 




Description 




set (A/b) 


: a'^b. Assigns a value to a variable (adds It to the 
library) 


+ 


ad(i(a,b) 


: a+t). Adds two VAUUEs 






: a-b- Subtracts two VALUEs 


* 




: a*b. Multiplies two VALUES a*b 


/ 


dlv(a,b) 


: a/b. Divides two VALUEs a/b bl=0 




delata <x) 


Removes a variable from the workspace 




sin (a) 


Calculates trigonomotric function sin of VALUE 




COS (a) 


Calculates trigonometric function cos of VALUE 




tan (a) 


Calculates trigonometric function tan of VALUE 




ctranapodd (X) 


Matrix transposed and complex conjugated 0 






Mabix transposed, (non-conjugated) (X.') 




tinea 


: X(i4)*Y(i j). Array multiply, element by element (.*) 






: t XOjyYO J)- Rislit array divide, element by element. (./ 
) 






: Y(I4)/X(ij). Left array divide, element by element. (A) 




power (3t,1Cj 


' X(l jj'^Xtl j). Array power, element by element (A) 




sqrt (X) 


Array square root, X(i j)>=0 




B\m (X) 


Sum of elements In vector x. 






: a= =b. Test - EquaK. TruQ=>1, fa]se=>0 




taatnec[(a/2}) 


: al =b. Test - NotEquaL. True=>l , fals6=>0 


> 




: a>b. Test If Higher. True=>1 » false=>0 






: a>=b. Test if Higher or Same, True^^l , false=>0 


< 


teaf:lo(a«b) 


: a<b. Test if Lower, True=>l , false=>o 






: a<sb. Teat if Lower or Same, f rue=>1 , falses>0 
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binary 


Function 


Description 


'-or I 


not <a) 


: -a. Logical NOT. 


& 




: a & b. Logical AND. 


I 
1 




• a A b Loaical OR. 






fiuiiri a veeJfir nf the araum^nts. 




getbuttXiixe (a) 


Gets the baseline no. a for the actual measurement 
from the database. Use It like this; bi-getbaseline(l) 




xnc3r (a) 


: a+1 , Increment a by 1 . 




g«tlxnea (xmiii. 


GetUnes returns a vector [Idx] of Indices to the desired 
frequer^Cles specified in vector [xd], for the (abclssa) 
vector specified by xmin,xmax,Nx. 

The function may b© used for frequency and time 
domain signatures. 




db20 (x) 


= 20*lQgiD(x/1). Returns the decibehvaiues of the 
elements in vector x, with reference to 1 (Unit). 




sigmoid (3C^ a, o) 


Sigmoid curve membership function. 

Returns a matrix which Is the sigmoid membership 
function evaluated at x (a canstant). pamms is a 2- 
element vector that determines the shape and position 
of this membership function. 

Specificaily, th© formula for this membership function is: 
slgmf(x, a, c) - 1-/(1 + exp(-a*(x-c))) 






Convert from Cartesian co-ordinates to polar co- 
ordinates. X and y are vectors. 

Returns a matrix where the first row is the r vector 
(fength), the second Is the f vector (angle). 

The Input vectors x and y express the abscissa and 
ondinate members. 
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binary 


7iinci2±on 


Description 




pol2cart (r,t) 


Suggested exchanging the existing 'poLcomplex' 
function. 

Convert from polar co-ordinates to Cartesian co- 
ordinated, 

Returns a matrix where the first row is the x vector 
(abscissa), the second is the y vector (ordinate). 

The Input vector r expresses the length scale, the 
vector t the angle Gcale, 




getmaaa (i) 


Suggested exchanging the existing 'getmeas* function 
without parameters, 

Gets the l^h measurement defined for the specific 
signature. 




g&tdi2iMii«as () 


Return the number of measurements assigned for this 
fault 






= x(lj). Get eJement in the l*th row and y column in 
matrix X. If x is a vector; may be left out 






= ZCx'")""". Tends against the highest values in vector x; 
more exclusively by increasing value of the constant 






= S(x^'^)"'"". Tends against the lowest values In vector x; 
more exclusively by increasing value of the constant m. 




£uz icygauB b (x r m , 
PCVB.X, val# axpo) 


Fuzzy weighting function. 

Ail input parameters are constants. 






Fundamental parameters of shaft orbit. 

Returns a vector with 2 elements where the first is the 
length of the major axis and the second is the fraction 
between the length of the m^or and the minor axis - 

The Input constants mx and phx express the magnitude 
and phase of the x-signal (horizontal) and my and phy 
express the magnitude and phase of the y-sfgnal 
(vertical). 
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binary 




Description 






Identification of shaft orbifs relation to an unbalance 
response. 

Express how large a portion of the estimated value at 
rpm that Is due to the quadratic term in the 
approximation, made by a second order polynomial tit. 

n/ec is vector of rpm values, cf is a vector of 
corresponding displacement values, and rpm is the 
speed where the frt \s tested. 






~ L x(|) . bum OT squares ot me eisniBnis in vector jc. 




vector trunc (Xf i / d ) 


Vector truncation. 

Returns a vector with the ith to Vnej'th elements from 
the vector x. 






p= max{x), Get8 the maximum value of the elements In 
the vector x. 




min (x) 


= mln{x). Gets the minimum value of the elements in 
the vector X. 




kurtosls (x) 


Calculate the kurtosis of the elements in the vector x, 




now 


Returns the currant date and time as a serial date 
number- FLOOR(NOW) is the current date and 
REM(N0W,1) is the current time, (Identical with the 
MATLAB command NOW) 






Returns the absolute value of the elements in the 

vector X, 






Returns the sign of the elements in the vector i.e. -1 
if the element is less than 0 else 1 . 




in«a&<x) 


= X//N. Returns the mean value of the elements In 
vector X. 






Retums the standard deviation of the elements in 
vector X. 




loglO (3C) 


= IoSid{x). Returns the logarithm with base 10 of the 
elements In vector x. 






t^Armraf<>4 VActnf i*if n noi'ni'<^ llnf^^riv Si:^9^fid between 

^ICI iDl Q,WS9 B VwwlUI " wuiiikw ill legally ^|-ica«^aw 

and including a and b* 






Returns the number of rows In matrix x. 




g6tclm<x# j) 


Returns the/th column from matrbc x. 



22 



binary 


]?uncti&n 


Description 




getrow(x,l) 


Returns the fth row from matrix x, 






Returns the number of columns \v\ matrix x. 




gatimaKvaluda (a) 


GetMaxValues returns the largest value of FUNCTION 
around tha abscissa values <L1NES> ± Dlines 




return (s/ as) 


Set the symptom & symptomstrength. 




tills t (ne&sno / 
dlxnno , a tart , flzid) 


Returns a vector with all elements on position efsmno 
In measurement meaeno from date and time sfert to 
date and time end. 




getBodA(l) 


Called from signature: OTAretum s 
v='[{l01]},{lOm.{<Oih{<02}JTo workspace: {RPMmiti, 
RPMmax, nRPM} ".from last run-up. 






Called from signature: GAP (scalar DC) 

Returns the co-ordinates of the Shaft Centre Line, [x,y]. 
Measurements are based on DC readings from two 
DPs mounted In arbitrary directions and the x,y co- 
ordinates are calculated using spatial transformation. 

To workspace; RPM 

Question; from which GAP sensor Is this callable and 
how is the other recognised?? 



